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Abstract 

We give a mathematical formulation for the Choi-Jamiolkowski 
(CJ) correspondence in the infinite-dimensional case close to one used 
in quantum information theory. We show that "unnormalized maxi- 
mally entangled state" and the corresponding analog of the Choi ma- 
trix can be defined rigorously as positive semidefinite forms on an 
appropriate dense subspace. The properties of these forms are dis- 
cussed in Sec. 2. In Sec. 3 we prove a version of a result from 
|14j characterizing the form corresponding to entanglement-breaking 
channel by giving precise definitions of the separable CJ form and 
the relevant integral. In Sec. 4,5 we obtain explicit expressions for CJ 
forms and operators defining Bosonic Gaussian channel. In particular, 
a condition for existence of the bounded CJ operator is given. 



1 Introduction 

In this paper we give a mathematical formulation for the Choi-Jamiolkowski 
(CJ) correspondence (|12j. [3]) in the infinite-dimensional case in the form 
close to one used in quantum information theory (see e.g. [13J). We show that 
there is no need to use a limiting procedure (cf. [6]) to define "unnormalized 
maximally entangled state" and the corresponding analog of the Choi matrix 
|3J since they can be defined rigorously as, in general, nonclosable forms on an 
appropriate dense subspace. The properties of these forms are discussed in 
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Sec. 2. An important question is: when the CJ form is given by a bounded 
operator. This is the case for entanglement-breaking channels: we prove 
this in Sec. 3 along with a version of a result from [H] characterizing CJ 
operators which correspond to such channels by giving precise definitions of 
a separable operator and a relevant integral. In Sec. 4 we obtain explicit 
expressions for CJ forms and operators defining a general Bosonic Gaussian 
channel. In Sec. 5 we give a decomposition of CJ form into product of the 
four principal types and a necessary and sufficient condition for existence of 
the bounded CJ operator. 

2 Positive semidefinite forms 

In what follows H, )C, . . . denote separable Hilbert spaces; denotes the 

Banach space of trace-class operators in TL, and &(T-L) - the convex subset 
of all density operators. We shall also call them states for brevity, having in 
mind that a density operator p uniquely determines a normal state on the 
algebra 25 (H) of all bounded operators in "H. Equipped with the trace- norm 
distance, ©("H) is a complete separable metric space. It is known [5], [I] 
that a sequence of quantum states {p n } converging to a state p in the weak 
operator topology converges to it in the trace norm. Moreover, it suffices 
that \im n (ip\p n \ip) = (ip\p\i/j) for ip in a dense linear subspace of H . 

Definition 1 A channel is a linear map ^:%('Ha) | — > %(T~Lb) with the prop- 
erties: 

1) $(&(Ha)) ^ &(Hb)] this implies that <E> is bounded map J^j and hence 
it is uniquely determined by the infinite matrix , where {\i)} is a 
fixed orthonormal basis in Ha- 

2) The block matrix [$ is positive semidefinite in the sense that 
for any finite collection of vectors {\ipi)} C % B 



The Choi-Jamiolkowski (CJ) form of the channel, associated with the 
basis {\i)} , is defined by the relation (0|) below. 

In 'Hb®'Ha we consider the dense domain which is invariant under "local" 
bounded operators Xb <8> Xa ' 

Hb x T-L A = hn {ipB <8>ipA'i>B e Kb, ?Pa e Ha} ■ 




(1) 
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Lemma 1 There is a unique sesquilinear positive semide finite form on 
Kb x Ha satisfying the relation 

^ (Vb®^;V>b®^) = (M®($a)($a\) Wb) = Wa\**(Wb)(M$a), 

(2) 

where \ip) = 

Proof. It is sufficient to show that for any \ipBA) = E Wb) ® Wa) e 

3 

E^i $ (i^)(^i)i^)>o, 

jk 

and \iPba) — implies that the above sum is equal to zero. Decomposing 

Wa) = ES%N)> we have \^ba) = E Wb) ® N) where I$b) = EciilV'e) 

* 3 

and the above sum is equal to J2{ip B \& Wb) whence the result follows. 

□ 

Note that for any orthonormal basis {(e^)} in ^ 

(e fc <8) V> A ; e fc ® ^) = {$aWa)- (3) 

The expression on the left is a natural form generalization of partial trace 
with respect to Hb- The relation fl3]) shows that for Q$ it is always given 
by the (form corresponding to) the identity operator I a- Similarly defined 
partial trace with respect to Ha not always exists; however this is the case 
when the expression $[Ia] is well-defined in the sense of [10] and then it is 
given by that operator. Positivity and the property ([3]) imply 

(lpB ® IpA, ^'b ® 1p' A ) I < [^<K (V>B ® ^B ® V'A) (V>B ® ^'a\ V>B ® 

< ll^flllll^HIWIII^II. (4) 

Conversely, any sesquilinear positive semidefinite form O on Kb x "H^ 
satisfying the condition (EJ) uniquely defines a channel $ such that fi$ = 
(via the reversed relation ([2])). The positivity of the form Vtq> can be used to 
prove the Kraus decomposition for $ similarly to the finite-dimensional case 
[T3] . Indeed, let Sj be the Hilbert space obtained by completion of Hb X Ha 
with respect to the inner product defined by the (factorized) form Then 
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any orthonormal basis in S) defines a countable collection of linear functionals 
on Kb x H A such that 

+00 

fi$ {tpB ® 0A! B ®^'a)=Ys h Wb ® ^) /i(^B®^). 

Z=l 

Define the linear operators VJ : "H^ — )■ "H_b by the relation ( , 0s|^z|V'a) = 
fi\4>B® 0a) j then 03]) implies that Vi are bounded operators satisfying 
EtT ^z*^ = j a and © implies the Kraus decomposition $ (p) = V^pVJ*. 

If the form fi$ is closable [15], then it is defined by the unique densely 
defined selfadjoint positive operator, which we also denote If the domain 
of the form is the whole He® Ha then the operator Q$ is bounded. The 
property (j3J) then reads 

Tr s O$ = I A . (5) 
In this case the defining relation ([2]) can be given a more familiar form 

Trfi $ (p <g> X) = Tr$ (p T ) X = Trp T $* (X) , p e %(H),X e*S(H), 

where T denotes transposition in the basis {\i)}, so that | "0a) ("0x1 = (I"0a) ("0a |) 
An example of nonclosable sesquilinear form is provided by the identity 
channel $ = Id, for which 

O w (0b ® V>a; 0b ® 0a) = (^b®V>a|0)Wb®^a), 
where (0| is the unbounded linear form on 7^ x Ha defined as 

+00 

<Q|-0i <8)-0 2 > = ^(#i>(^2>; Vi,^2 e ^a, 
i=i 

and is the dual antilinear form which represents "unnormalized maxi- 
mally entangled state", \Q) = Xl^T K) ® K)- The relation 

tt$ = ($ <g> Id A ) (fi w ) 

holds in the weak sense i.e. as equality for the forms defined on Hb x Ha- 
Notice that (X4 <g> X B ) \Q) = ( I ® X B XT) If is another basis 

such that = £7|i) for a unitary [/, then = ( I ® /7/7 T ) and 

0' $ = (J ® /7/7 T ) (J ® /7/7 T )* . 
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In the case of bounded f2$ this implies that is the same for all choices 

of the basis {\i)} ■ Notice also that transposition in the basis {\i)'} is given 
byX T ' = (UU T ) X T (UU T )* . 

Fix a state a in &(Ha) of full rank, and let De the basis of eigen- 

vectors of a with the corresponding (positive) eigenvalues {Xi}^. Consider 
the purifying vector 

+00 

= X> 1/2 io ® 10 

i=l 

in the space %a ® Ha- Then the state 

pz(a) = ($®kU)(|Vv)(?/v|) e G(H B ®H A ) (6) 
satisfying TtbP<s>{<j) = <y uniquely determines the channel $ via the relation 

= A- 1/2 A- 1/2 Tr A (/ B ® b'XilWCo-), 
see [TT]. The connection between p$(cx) and 0$ is 

(iI) B ®iI)a\p${p)Wb®^a) = fi * ® v 1/2 *Pa; $ b ® v 1/2 *P'a) ■ 

Note that 0$ is uniquely defined by its values on the dense domain T> = 
H B x (T l/2 (U A ) due to the property flU). 



3 Separable operators and entanglement-breaking 
channels 

Let a be a state in &(Ha) of full rank and Q is a bounded positive operator 
satisfying (jSJ), then ct^a = ® cr 1 ^ 2 ) [Ib <8> cr 1//2 ) is a density operator in 
T-Lb^Wa such that Trg erg a = a. Let us remind that a state p G &(Hb®'Ha) 
is called separable if it is in the convex closure (in the weak operator topology 
and hence in the trace norm) of the set of all product states. Separable states 
are precisely those which admit the representation 

p= (p B (x) ® p A (x)) /i(dx), (7) 
Jx 

where p, is a Borel probability measure on a complete separable metric space 
X, see PI]. 
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A channel $ is called entanglement-breaking if for arbitrary state u> G 
& (Ha® Ha) the state ($Cg>IcU)(w) is separable. Channel $ is entanglement- 
breaking if and only if there is a complete separable metric space X, a Borel 
&(Hb) -valued function x !->■ Pb(x) and a probability operator- valued Borel 
measure (POVM) M A (dx) on in Ha such that 

$(P) = J p B (x)p p (dx), (8) 
.V 

where p p (E) = Tr pM A (E) for all Borel i£ C X . If p = p$(cr), then Pb(x) in 
([8]) is the same as in ([7]) while M^cte) is defined by the relation 

(i// a \M a {E)\*I>a) = [ {v- l/2 ^ A \p~Zx)\c?- l / 2 4, A )p{dx)- ip A ^' A G a 1 ' 2 (Ha) , 

J E 

where the complex conjugate is in the basis 

Definition 2 A bounded positive operator Q in Hb <8> H A satisfying is 
called separable if it belongs to the closure, in the weak operator topology, of 
the convex set of operators of the form Pa ® M a , where {M a } is a finite 

a 

resolution of the identity in Ha and p a G &(Hb)- The weakly closed convex 
set of separable operators will be denoted CCba- 

Lemma 2 For Q G £ba the operator norm < 1. 

Proof. It is sufficient to prove it for Q = J^Pa <8> M a , since the weak 

a 

operator limit does not increase the norm. Then 

||fi|| = sup {il>\£l\ip) < sup V/b ® M a \ip) < (ip\ip) = I. 

11-011 = 1 11-011=1 a 

□ 

It follows that in the definition [2] it is sufficient to consider sequences of 
operators, weakly convergent on a dense subspace of %b ® Ha- 

Equipped the definition [2] and the construction of the integral (GO) below 
we can prove a rigorous version the corresponding result from [T3] . 
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Proposition 1 If the channel $ is entanglement-breaking then its form is 
given by a bounded operator Q$> G £ba-IJ & has representation then 



r2$ 



/ 



p B (x) ®M A (dx), 



(9) 



x 



where the integral is defined in the proof. 

Conversely, if a bounded operator Q G £ba then Q = where the 
channel $ is entanglement-breaking. 

Proof. The proof of the first statement requires some theory of integration 
with respect to a POVM. 

Let X be a complete separable metric space with a— algebra of Borel sub- 
sets £>, let % is a separable Hilbert space and {M(E); B G B} a POVM on X 
with values in 23 ("H). Then for any ip G "H the set function {(i[)\M(E)\i/;); B G B} 
is positive finite measure with total variation H^H 2 ; and for any G H 
the set function {(i[)\M(E)\i/;'); B G B} is a complex measure of finite total 
variation \\if)\\ as follows from the inequality 



due to positivity of the operator M(E) for arbitrary Borel E C X . 

A function x — > p(x) with values in ©("H) will be called Borel function 
if the scalar functions x — > are Borel for all G "H. Let 

p be a a— finite measure on X then the function x —> p(x) will be called 
measurable if the functions x — > (^|p(:c)| ■?/>') are //—measurable for all ■?/>, ip' G 
"H. This implies that for any A G 23 ("H) the scalar function x — > Trp(x)A 
is /i— measurable. Since &{H) is separable with respect to the trace norm 
distance, this implies, by theorem 3.5.3 from [7], that x — > p(x) is strongly 
measurable in the sense that there is a sequence of simple Borel functions 
x — > p n (x) such that \\p(x) — p n (%)\\i (mod/i). If p is a probability 
measure then the Bochner integral J x p(x)p(dx) = \im n ^ ^ J x p n (x)p(dx) 
exists and is an element of (5("H). 

Now let $ be entanglement-breaking channel with the representation ([8]). 
If x — > Pb{x) is a Borel function with values in &{T-Lb) and M^(<ix) a POVM 
in Ha we define the integral J ps(x)®MA(dx) G *B('Hb®'Ha) as follows. For 



a simple function p n (x) = Pa^E a (x), where {E a } is a finite decomposition 



\(i>\M(E)W)\ < l - [c(m(EM + c-'imiEM')} ; 



c> 0, 



(10) 



A' 



a 
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of X , we put Q n = f p n (x) g) M A (dx) = J2p a ® Ma(£? ). Apparently 
f2 n G £ba, hence ||f) n || < 1. We also have 

(lp B ®1pA\ttn\lpB®'>pA) = j \^B\pn{x)W B ){^ A WA{dx)\i> A ), (11) 

X 

for all ipBj iPb ^Hb, ^a^'a e Ha- For two simple functions p n (x), p' m (x) we 
have by the inequality (jTOl) 

< 5 N>b II IIV-bII y \\Pn{x) - p' m {x)\\ x [c$ A \M A {dx)\$ A )+ C -\$ A \M A {dx)W A )] ■ 

X 

Take {p n (x)} such that ||p(x) — p n (a;)|| 1 — >■ (mod/v) where a is a state in 

&(H A ) of full rank, then it follows that the forms (ipsA I^IV'ba) converge 

on a dense subspace of %b ® Ha and are uniformly bounded. Hence there 

exists uniquely defined bounded operator Q, which is the limit of Q n in the 

weak operator topology. We define f Pb(x) <S) M A {dx) = Q, and it follows 

x 

from the definition that Q G £ba- 

From the definition ([2j) of the form Q$ we obtain 

n*(V'B®Vu;V'B®^A> = J ^B\pB{x)W B ){4)' A \M A {dx)\$ A ), 

X 

which is the limit of fill I) . Since, on the other hand, (TTT1 ) converge to the 
form defined by the operator Q, we obtain that the form Q<p is defined by 
the operator Q. 

Conversely, let a bounded operator fl G £ B a- Fix a state a in &(H A ) of 
full rank. Using the fact that weak convergence of operators Q = ^ p a <S> 

M Q G implies weak convergence of density operators oba = Yl Pa ® 

a 1/2 M Q o- 1/2 , we can prove that the state cxba = (/# ® a 1/2 ) (Jb <8> o" 1/2 ) is 
separable. Basing on the decomposition ([7j) for cr and arguing as in [11] we 
can construct POVM M A (dx) and the family of states Pb(x) such that f2 is 
given by (Q and hence Q = Q$ for the corresponding entanglement-breaking 
channel $ given by (jSJ). □ 
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4 Bosonic Gaussian channels 



Let (Z, A) be a coordinate symplectic space, dim Z = 2s, with the symplectic 
form 

A(z, z ) = z f Az'; A = diag J ~ 

Jj=l,...,s 

where * denotes transposition in Z. Let % be the space of irreducible repre- 
sentation of the Canonical Commutation Relations (CCR) 

W(z)W(z') = exp (^ A ( z > z ')) W ( z + z ')- ( 12 ) 

Here W(z) = exp(iR-z); z G Z, is the Weyl system and R = [q\,pi, ■ ■ ■ , q s ,Ps] 
is the row- vector of the canonical variables in "H, see [8] for more detail. We 
will continue to denote T transposition in "H associated to a basis from 
(TT2"|) it follows that W(z) T = exp(iR T ■ z); z G Z, is the Weyl system for 
the symplectic space [Z, — A). The canonical transposition associated with 
the Fock basis in % is given by R T = [qi, —p\\ . . . , q s , —p s )] all the others are 
obtained from this by unitary conjugations. In what follows transposition is 
arbitrary if not stated otherwise. 

Let Ha be the representation space of CCR with a coordinate symplectic 
space {Za-i A^) and the Weyl system Wa{-), with a similar description for 
H,b- Let $ : T("Hyi) i — y T('Hb) be a channel. Then the following relation 
holds 

ft* = j^— [ W B (-z) ® $* (^ B (z)) T d 2s z, (13) 

in the sense of forms defined on Hb x Ha- Here d 2s z is the element of sym- 
plectic volume in Zb and T is a transposition in T-La- 
Indeed, 

(W B (z)) = M(|$$'|)Wb(s), (14) 

where is trace-class operator and hence (fT4|) is square-integrable 

as a function of z G Zg. Similarly, the function (0|W(— z)^') is also square- 
integrable. By the inversion formula for the noncommutative Fourier trans- 
form (cf. ch. V of [5]), applied to the right-hand side of (THj) . 

ft* {1>b ® *Pa; ® = (V'b|$(|^>(^|)|V'b) 

l 



(27T) 



<V>s|Wb(-«)|^b)(^I^ (Wb(s)) I^a)^ 
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so that finally we get the formula ( I13|) . 

Now let $ be a (centered) Gaussian channel [1], 

$* (W B (z)) = W A {Kz) exp (-^z'pz) 



zez 



B, 



where 



H > ± l -A K ; A K = A B - K l A A K. (15) 



Then ( TT3|) implies 

1 /' T / I 



n * = j^ky J Wb (-*) ® w *( Kz ) 1 ex p I ~i^H d2Sz - (16) 

The unitary operators Wba(z) = W B (z) <g> Wa(—Kz) t = exp(iRsA • z) , 
where 

Rba = R b ®I -I® RaK, 

satisfy the Weyl-Segal CCR with (possibly degenerate) symplectic form de- 
termined by the matrix Ak- This representation in the space 7i B <S> Ha is 
reducible (even if Ak is nondegenerate, see the footnote below). Finally 

fi* = j^y J exp {—iRba ■ z) exp \^\z l \iz^ d 2s z. (17) 

If ji is nondegenerate, then the form is given by bounded operator 
with 

< — 1=. (18) 
Vdet \i 

This follows from ( ITT)) taking into account the fact that norm of the Weyl 
operators is equal to 1. 

In what follows we shall consider the four basic cases. Later it will be 
convenient to use reverse enumeration, so we start with the last, the most 
degenerate case. 

Case 4. If fi = then Ak = and R B a is the vector operator with 
commuting selfadjoint components. Thus 

= TTTTI / ex P (~ iR BA ■ z) d 2s z = 5 (R BA ) , (19) 
(2n) s J 
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where 8 (•) is Dirac's delta-function. In this case Q$ is a nonclosable form. 
Note that for the ideal channel (A = B) this gives 

\Q)(Q\ = J— [ W B (-z)®W A (z) T d 2s z = 5(R BA ), (20) 

where R BA = Rb ® I - I ® 

Case 3. If /x > 0, A K = then again R BA = R B <g> / - 7 <g> -R^K" is 
the vector operator with commuting selfadjoint components and the integral 
( I17p is just the multivariate Gaussian density as a function of R BA '■ 

Since the spectrum of R BA contains 0, we have = ^ ct ■ 

Proposition 2 Assume that Ak is nondegenerate, then 

\n^\ = -== 1 (21) 

Met A K det [abs (A^V) + 1/2] 

Case 1. If jj,— \ Ak is nondegenerate, then 

fi* = , 1 = exp ( R BA tR BA j , (22) 
y^det (/i-iAj 

where e = arccot (2A^- 1 /i) A^ 1 . 

Case 2. If ji — |A# is maximally degenerate i.e. rank (/i — |A^J = s, 
then = ^= Ak ^ ' w ^ ere Po is the projection onto the kernel of positive 
selfadjoint operator R BA fi~ 1 R t BA — si and ||f2$|| = ^ dc * a a . ■ 

Proof. If Ajc is nondegenerate then \i is also nondegenerate, so is 
given by bounded operator. Rewriting ( Tl6l) and again using the inversion 
formula we get 

^ = (2 7 r)VdetA / , / 6XP { ~ IRbA • Z) 6XP (4^) ^< Z = 7^ PK ' 

where d\ s K z = y/det Axd 2s z is the volume element corresponding to the sym- 
plectic form z t Ajcz', and px has the expression in the canonical variables R BA 
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as the Gaussian density operator^ with zero mean and the covariance matrix 
\i. The value ( 12"T|) is just the maximal eigenvalue of this Gaussian density op- 
erator, multiplied by (det A^) ' . There is a nondegenerate transformation 
T such that 



fx = T l pT = diag 



Hj 
/Mj 



A = T'A X T = diag 



-1 

1 



(23) 

where fij > |, j — 1, . . . , s (see the next Section). Then /\r l p = diag 

and A^/i = TA~ 1 /iT~ 1 is matrix of the operator with eigenvalues ±ipj, so 

^ T _1 has the eigenvalues u 7 - 
Hj J 

of multiplicity 2 . The operator px splits into the normal modes decomposi- 
tion 



that the operator abs (Ajf/x) = Tdiag 



Pk 



(24) 



with p^> being the elementary one-mode Gaussian density operator 



P 



U) 



Pj +2 V^j + 2 



2 



(25) 



where = | (g| +Pj — l) is the number operator for the j— th mode (see 
ch. V of [8]). Here the new canonical variables R = [qi, ■ ■ ■ ,p a ] are related 
to the old ones by the formula R = RbaT. The maximal eigenvalue of px is 
thus equal to 



5* + 



det [abs (A^V) + V 2 ] 



Since /x — |A^ = A^ (A^/i — |J , the condition that // — |A# is nonde- 
generate (Case 1) is equivalent to p,j > |, j = 1, . . . , s, i.e. the decomposi- 
tion ( 12~4"|) has no pure component. Coming back from (1241) . (|25p to the initial 



1 Notice that is not a proper density operator in the space ® Ha since Rba 
generate a reducible representation Wba(z) of CCR in that space. Actually pk is tensor 
product of the Gaussian density operator in the space where Wba{z) act irreducibly with 
the identity in the complementary space, reflecting the multiplicity of the representation. 



12 



canonical observables Rba gives 

p K = cexp (-RbacRba) > ( 26 ) 

where 

c = Il^^= . 1 (27) 



and e is found from 

2A x V = coteA^, (28) 

whence the formula (|22j) . 

The case 2 corresponds to pj = |, j = 1, . . . , s. Then px is the projection 
onto the kernel of the positive selfadjoint operator 

s s 
3=1 3=1 

so that 

n * = / a l a g o (RbaP^Rba - si) 
Vdet A^ 

(Kronecker's delta), and ||fi$|| = ^ de ^ A ^ .D 

Example. Consider attenuator /amplifier in one mode, 



A 

given by 



-1 

1 



/ 777 

<F (W(z)) = W{kz) exp |z| 



z G M 2 . 



where k,m > 0. Then ffl5|) reduces to m > - — 2 — an d det A^ = (fc 2 — 1) 
The relation (T2~Tj) gives 



in 



m + 1 2 1 



The channel is entanglement-breaking if and only if m > i.e. m + 

|fc 2 -l| l/c 2 — 1| 

- — 2 — - > max{l,/c 2 } [H] which agrees with Lemma 1. Also if 1 — - — ^ — - < 
m < then the channel is not entanglement-breaking while still < 1. 
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\k 2 -l\ 

Together with m > - — 5 — this gives k > 1 (amplifier) and the lower bound 



m > max 



3-fc 2 k 2 - 
2 > 2 

T 



Assuming the canonical transposition we have = [q; —p]. Then, 

( lo- 
using the relation (1231) . one can obtain in the case 1 m > - — =- 



1 , m + 



k 2 -i 



n * = I 6XP 1 " 2IF-1I ln ~ + (PB + kpA 



(29) 

|fc 2 -i| _i 
The case 2 corresponds to m = J — 5— then one obtains = |A; 2 — 1| P > 



where Po is the projection onto the eigenspace of the operator (qs — kqA) + 
(ps + kpA) 2 , corresponding to its lowest eigenvalue \k 2 — 1| . 

In general, the decomposition ( 1241) means that in the case 1 the operator 
fi$ can be decomposed into tensor product of operators of the form (1291) . and 
similarly in the case 2. 



5 A decomposition of the Gaussian CJ form 



Recall that 2s = dim Zb and denote by r a = rank a - the rank of a 2s x 
2s— matrix a. The following result is a generalization of the Williamson's 
lemma (cf. [2]). 

Lemma 3 Let fx be a real symmetric matrix, Ak - a real skew-symmetric 
matrix such that fx — \Ak > 0. Then there is a nondegenerate matrix T such 
that 



where 



diag 



T l A K T 



-1 

1 



fx 

1/2 



A 







} r A K 

K - r A A - 
}2s - r„ 



(30) 
(31) 



jx = diag 



j=l,...,r AK /2 



fij 
fij 



i=l,...,r A ,,/2 
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and fij > 1/2. 



Notice that r„ can be odd. Denote 



di 



2s 



r M the 



dimensionalities of the last two blocks in the decompositions f l30|) . (l3Tj) . Let 
us further arrange the block diagonal matrix jl as 







}d 2 



by putting first the blocks with /ij > 1/2 and then - the blocks with fij = 1/2. 



We have r^ K = di + d 2 ,r 



^1 + ^2/2 + ^3, whence 



= r AR - - 2(r, 



d 2 = 2(r, 



Let ej = T~ 1 ej]j = 1, . . . , 2s be the basis in Zb in which /i, A^- have the block 
diagonal form fl30|) . fl3~TT) and let Z fc be the d k — dimensional subspace spanned 
the vectors Cj corresponding to the k— th block in the decompositions, k = 
1, . . . , 4. Then we have the direct sum decomposition 



Zi + Z 2 + Z3 + Z4 



(32) 



By making the substitution T l z = z in ( fl7|) . we have z = [zi,z 2 ,z 3 , z^\ 
and 

= (2tt) s |detT| / / / / 6XP ^ (- iR BATz k - ~4/* (fe) %^ dz x dz 2 dz z dzL 

where p,^) — ^d 3 /2, A(4) = 0d 4 and the components of RbaTz^ and RbaTzi 
commute for k ^ I by fl3Tl) . Hence the exponent under the integral splits 
into product of four mutually commuting exponents, and the CJ form 
can be decomposed into the product of commuting expressions of the types 
considered in the cases 1-4 above (with possibly odd dimensionalities for 
zi,z 2 ): 

^$ = Yl J. ex P y-iRBATh - -z{fi{k)ZkJ dz k . (33) 
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This product can be further transformed into tensor product in the space 
%b ® Ha as follows. Consider the direct sum Z B + Z A equipped with the 
symplectic form defined by the skew-symmetric matrix 



A 



AB 



' I 


K* ' 




' A B 







" I 


' 




A K 


-K*A A 





I 







-A A _ 




K 


/ 




_ ~A A K 





This form is nondegenerate since it is determined by the product of nonde- 
generate matrices; moreover, its restriction to Z B coincides with Ak- The 
unitary operators 



W(z B , z A ) = exp i (R B az b - I © R A z A ) = exp % I [R B , -R 



I 
K I 



z B 

ZA 



give an irreducible representation of CCR on the symplectic space (Z B + Z A , A AB ) 
and at the same time W(z B ,0) = exp i (R BA z B ) = W B (z B ). It follows that 
the direct sum decomposition (|32|) can be extended to the decomposition 



Z B + Zj 



z' 3 ®z' 4 



Zq, 



(34) 



where Z' 3 D Z 3 , Z' 4 D Z 4 and © means that A AB is nondegenerate on each 
of the five subspaces (provided they are nontrivial) and zero between the 
different subspaces. 

To prove this, we use the fact: if (Z, A) is a symplectic space and Z a 
subspace of Z such that the restriction A\% is nondegenerate, then there is 
a unique subspace Z 1 - such that Z = Z © Z 1 - . From the decomposition ([311) 
and the fact that A AB \ Zg = A K it then follows that Z B + Z A = Z\ © Z 2 © M, 
where M D Z3+Z4. Then (M, A ab \m) is itself a symplectic space and A ab \m 
is identically zero on Z3 + Z±. The basis {e-,} in Z% + Z4 can be comple- 
mented by the system {hj} in M such that A AB (ej, hk) = Sjk, A AB (ej, e^) = 
A AB (hj, hk) = 0. Let the subspaces Z' 3l Z' A C M be spanned by the vec- 
tors ej,hj such that the corresponding ej span, respectively, Z 3) Z A . Then 
by construction A AB is nondegenerate on Z' 3 , Z 4 and equals zero between 



Zx, Z 2 , ZL Zi. Denoting Z 



Z' © Z' 



we obtain 



Z\ © Z2 ^3 ^4 

But this decomposition means that H, B © T-L A can be splits into tensor 
product T-Li®%2® H3 © %4 © Ho such that exp (iR BA Tzk) acts nontrivially 
in %k for k — 1, . . . , 4 so that the product fl33|) can be transformed into the 
tensor product. 

Since in the cases 1-3 the integrals in the product are given by bounded 
operators, we can complement (|18p as follows: 
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Proposition 3 Nondegeneracy of the matrix /i is necessary and sufficient 
for the form f2$ to be defined by a bounded operator. 

Finally, let us give interpretation of the decomposition (132]) in terms of 
open system dynamics. Then we have the composite system Za®Zd = Zb® 
Ze (input+noise (distortion)=output+environment) which evolves reversibly 
according to the unitary operator U in the space Ha ®Hd = He® He- The 
dynamical equations of the channel in the Heisenberg picture can be written 
as 

R' B = U* (R B ®I e )U = R a K ®I d + I®R d , (35) 

where Rd is the vector of noise variables having the commutator matrix 
Aj<- = Ab — K*AaK and the covariance matrix \x. The lemma [3] implies that 
Rd = [Rdi RdV^' 1 where R q D is the d± + d 2 = r& K — dimensional subvector of 
quantum noise canonical observables with the commutator matrix A and the 
covariance matrix jx and R C D is d% + d± = 2s — ta k — dimensional subvector 
of commuting classical noise variables (which commute also with R%). The 
summands in the decomposition (|32l) correspond to 1) quantum noise ob- 
servables in the nondegenerate Gaussian state, 2) quantum noise observables 
in the pure Gaussian state, 3) classical noise variables with positive variance, 
4) trivial classical noise variables with zero variance. 

Acknowledgments. The author is grateful to V. Giovannetti and M. 
E. Shirokov for discussions. The first version of this work was partially 
supported by RFBR grant 09-01-00424 and the RAS program "Mathematical 
control theory" . The new improved and extended version was written during 
the author's stay at the Mittag-Leffier Institute in the frame of the program 
"Quantum Information Theory", Fall 2010. 

References 

[1] Caruso F., Eisert J., Giovannetti V., Holevo A. S., Multi-mode 
bosonic Gaussian channels, New Journal of Physics 10, 083030 2008; 
larXiv:0804.05TT[ quant-ph] . 

[2] Caruso F., Eisert J., Giovannetti V., Holevo A. S., The optimal unitary 
dilation for bosonic Gaussian channels; larXiv: 1009. 110 8 [quant-ph] 



17 



[3] Choi M.-D., Completely positive maps on complex matrices, Linear Alg. 
and Its Appl. 10, 285-290, 1975. 

[4] Davies E. B., Quantum theory of open systems, Academic Press, London 
1976. 

[5] Dell' Antonio G.F., On the limits of sequences of states, Commun. Pure 
Appl. Math. 20, 413-430, 1967. 

[6] Giedke C, Cirac J. C, The characterization of Gaussian operations 
and Distillation of Gaussian States , Phys. Rev. A 66, 032316, 2002; 
quant-ph/0204085| 



[7] Hille E., Phillips R., Functional analysis and semigroups IL 1962. 

Holevo A.S., Probabilistic and statistical aspects of quantum theory, 
North Holland, Amsterdam 1982. 

[9] Holevo A.S., Entanglement-breaking channels in infinite dimen- 
sions, Problems of Information Transmission, 44, n. 3, 3-18, 2008. 
larXiv:0802.0235{ quant-ph] . 

[10] Holevo A. S., The entropy gain of infinite-dimensional quantum chan- 
nels, larXiv: 1003.576 5 [math-ph] . 

[11] Holevo A. S., Shirokov M. E., Werner R. F., Separability and 
entanglement-breaking in infinite dimensions, Russian Math. Surveys, 
60, N2, 2005; |arXiv:quant-ph/0504204"l 

[12] Jamiolkowski A., Linear transformations which preserve trace and posi- 
tive semidefiniteness of operators, Reports on Mathematical Physics, 3, 
N4, 275-278, 1972. 

[13] Nielsen M. A., Chuang I., Quantum Computation and Quantum Infor- 
mation, Cambridge University Press 2000. 

[14] Owari M., Plenio M. B., Polzik E. S., Serafini A., Wolf M. M., 
Squeezing the limit: quantum benchmarks for the teleportation 
and storage of squeezed states. New J. Phys. 10, 113014, (2008); 
larXiv:0808.2 260 [quant-ph] . 



18 



[15] Reed M., Simon B., Methods of Modern Mathematical Physics. Vol.1. 
Functional Analysis. Academic Press 1980. 



19 



